Resonance energy transfer near metal nanostructures mediated by surface plasmons 
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We develop a unified theory of plasmon-assisted resonance energy transfer (RET) between 
molecules near a metal nanostructure that maintains energy balance between transfer, dissipation, 
and radiation. We show that in a wide range of parameters, including in the near field, RET is 
dominated by plasmon- enhanced radiative transfer (PERT) rather than by a nonradiative transfer 
mechanism. Our numerical calculations performed for molecules near the Ag nanoparticle indicate 
that RET magnitude is highly sensitive to molecules' positions. 



I. INTRODUCTION 

Resonance energy transfer (RET) between spatially 
separated molecules 1 ^ plays an important role in diverse 
phenomena across physics, chemistry and biology. Exam- 
ples include photosynthesis, exciton transfer in molecu- 
lar aggregates, energy exchange between proteins,— i 4 - and, 
more recently, between excitons in quantum dots (QDs) 5 
and in QD-protein assemblies.-^ During past decade, sig- 
nificant advances were made in RET enhancement and 
control by placing molecules or QDs in microcavities^— 
or near metal films and nanoparticles (NPs)J£~— The 
coupling between molecular dipoles and surface plasmons 
(SP) in metal opens up new RET channels. The ability to 
control RET rates by adjusting dipoles' positions relative 
to metal surface is important in biomedical applications 19 
such as, e.g., SP biosensors^ 

Near plasmonic nanostructure, RET from a donor to 
an acceptor is governed by the interplay between several 
processes. The energy of the excited donor can either 
be radiated, dissipated, or absorbed by the acceptor and 
each of these channels is affected by the nearby metal in 
its own way. In a closely related phenomenon, plasmon- 
enhanced fluorescence, the decay rates in nonradiative 
and radiative channels depend differently on the distance 
between molecule and metal surface, d. The measured 
fluorescence^ 5 . - — from molecules attached to a metal NP 
indeed shows that with decreasing d, SP enhancement is 
followed by quenching, in agreement with theoryi2i~— A 
similar, albeit somewhat more complicated, scenario is 
expected when a donor and an acceptor molecules are 
placed nearby a plasmonic nanostructure; i.e., the en- 
ergy transfer from the donor to the acceptor should be 
strongly affected by dissipation in metal and by plasmon- 
enhanced radiation. However, no RET theory including 
all relevant energy flow channels has yet been available. 
It is our goal to provide such a theory here. 

To highlight the issue, recall the famous Forster's 
formula for energy W^ d transferred from donor to 
acceptor— ^ 

where W d is the donor's radiated energy, is its spec- 



tral function, o~ a (uj) is the acceptor's absorption cross- 
section, D ad is the dipoles' electromagnetic coupling at 
distance r ad and k is the wavevector of light. In the 
near field (kr ad < I), we have D ad = q a d/rl d {q a d is the 
orientational factor) and RET changes with distance as 
(rp/fad) , where rp is Forster's radius. In the far field 
(kr ad ^> 1), RET is dominated by radiative coupling 
l-^adl ^ k 2 / r ad leading to weaker r~ d dependences^ 
Eq. ([!]) is derived from first-order transition probability 
under the perturbation D® d . 

For molecules near a plasmonic nanostructure, Eq. ([T]) 
must be modified. The standard model by Gersten and 
Nitza n 30 ! 31 and its extensions to planar and composite 
systems 3 ^— incorporate SP in the transition's interme- 
diate states, and thus Eq. ([TJ still holds albeit with new 
coupling D ad which now includes SP-mediated channels. 
However, this model accounts for neither dissipation in 
metal nor plasmon-enhanced radiation channels and, as 
a result, yields enormous (up to I0 5 ) RET enhancement 
that contrasts sharply with the much more modest (~ 10) 
increase ^ 14 ' 16 " — and even reductio n 10 ' 15 of measured 
RET rates. 

Here we present a unified theory for RET near metal 
nanostructures based on the classical approach that ac- 
counts accurately for the full energy flow in the system. 
We show that Eq. (1) is replaced with 

W ad 9 f dw H - - 2 

where j d is the donor's free space radiative decay rate, Td 
is its full decay rate, fd and a a are the modified spectral 
function and absorption cross section, respectively, and 
the coupling D da includes high-order SP-assisted tran- 
sitions. For a low-yield donor, j d should be replaced 
with the free space fluorescence rate j d . We also iden- 
tify plasmon-enhanced radiative transfer (PERT) as the 
dominant RET mechanism in a wide parameter range. In 
the far field, we extract from Eq. ((2|) a general formula 
for PERT from remote donors to an acceptor near the 
metal surface that extends radiative RET theory 4 ^ to 
plasmonic systems. In the near field, our numerical calcu- 
lations of RET near Ag NP (see inset in Fig. [IJ show that 
PERT is the dominant mechanism here as well. Depend- 
ing on system geometry, RET can either be enhanced 
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FIG. 1: (Color online) RET vs distance for R = 30 nm Ag NP 
is shown at 9 — it/3 with (a) d a = dd — d and (b) d a = 20 nm, 
dd = d using the full Eq. @, the nonradiative (NR) channel 
only, Forster's transfer Eq. {l}, and the Gersten-Nitzan (GN) 
modelM 



or reduced as compared to Forster's transfer, consistent 
with experiment lo- 



tion with dipole sources^ 

AlTLU 2 



E(r,w) 
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where G(r, r';u>) is Maxwell's equation Green's dyadic, 
satisfying V x V x G — e(r, w)(w/c) 2 G = I, and e(r, w) 
equals metal permittivity, e(u>), inside the metal region, 
and that of the outside medium, eo, otherwise. The quan- 
tity of interest is energy absorbed by the acceptor in the 
unit frequency interval, 



= 1111 [Pa( W ) • E ( r a, W) = - 
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where we used E(r a ,w) = p a (u)/a a (ui) from Eq. A 
closed system for Pj(uj) is obtained by using Eq. to 
eliminate the electric field from Eq. ([3]), 



Pj(uj) + otj ^2 D jk{u)Pk{u) = S jd p° d (uj), 

k 

where we introduce the frequency-dependent matrix 
Aituj 2 



D jk {u>) = 
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Expressing p a from Eq. ©, we obtain 
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where D a d = D a d [1 - adDd a a a D a d] is donor-acceptor 
coupling that includes high-order transitions, and 



1 + D 3j {u)a 3 {u) 



(9) 



is the molecule's dressed polarizability satisfying the re- 
lation 



II. THEORY OF PLASMON-ASSISTED 
RESONANCE ENERGY TRANSFER 



We consider a donor and an acceptor near the sur- 
face of a metal nanostructure which are represented by 
pointlike dipoles located at Vj with induced moments 
Pj(u>) = pj(ui)ej oriented along ej (j — a,d). The dipoles 
are driven by the common electric field, 



Pj(u)) = a j (w)E(r j , u) + 5 jd p° d (u), 



(3) 



where a.j(oj) = a'j(oj) + ia"(oj) is complex polarizabil- 
ity assumed here isotropic, p d (w) = a d (u})e d Eo is the 
donor's initial dipole moment with some constant Eo de- 
pending on excitation, and 5j k is Kroniker's symbol. The 
electric field E is, in turn, the solution of Maxwell's equa- 



(10) 



which expresses the energy balance between total extinc- 
tion described by a", external losses such as radiation 
and dissipation in metal encoded in D"j(w), and absorp- 
tion in the presence of environment (right-hand side). 

To gain more insight, recover first Forster's RET from 
Eq. flHJl. For a high-yield donor (a d = 0), Eq. ([10]) yields 



the optical theorem a dQ 
a j0 = 



\k 3 \ 
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2 , where 
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(11) 



is polarizability in radiation field and we use free space 
expression for = — i^k 3 . The near field coupling is 
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with r = r/r, while a a D® a ~ a a k 3 is negligible. The 
radiated energy of an isolated donor can be derived in a 
similar manner as 



W d 



El 



(13) 



Using the optical theorem, Eq. © leads to Eq. ([T]) with 

= ^*a>). M")= rf^n) >v (14) 
3 J duW d ' (w) 



where the free space donor's spectral function /d(w) is 
integral-normalized to unity. 

Turning to the general case, we note that for a high- 
yield donor, the energy balance relation Eq. (fTUf implies 
the optical theorem in an absorptive environment, 
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where Fj = —fj^D'Jj is the molecule's /mH decay rate^ 
and Yj — |^ 3 M^ i s its radiative decay rate (fij is the 
dipole matrix element). Using this relation and normal- 
izing Eq. (jHJ to the radiated energy of an isolated donor 
[Eq. ([15])]. we obtain 



1 dW n 
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W d 87rfc 4 r d (w) 
which leads to Eq. ([2]) after frequency integration. Here 



Airk 



(17) 



are the acceptor's absorption cross section and the 
donor's spectral function modified by the environment 
[compare to Eq. fH|) ]. Note that, in the presence of 
metal, /d(w) is no longer integral-normalized to unity. 

Equation (I16[) includes all relevant energy flow chan- 
nels in the system. Interactions of the molecules with the 
metal alter the positions and shapes of the optical bands. 
While the coupling D a d is enhanced due to plasmon- 
mediated channels, the factor 7d/Td accounts for RET 
quenching due to the donor's energy transfer to the metal 
followed by dissipation and radiation. The absence of 
this factor leads to spuriously large RET^2r— Note that 
Eq. ([2]) was obtained for a high-yield donor with no as- 
sumptions on molecules' emission or absorption spectral 
bands, which are usually broad and asymmetric due to 
vibrational and rotational modes. Rigorous treatment 
of molecules' internal relaxation processes would require 
fully quantum-mechanical consideration which is beyond 
our scope. However, if we assume Lorenzian lineshape 
for the donor's effective polarizability a^u;), which is a 
reasonable approximation in most cases, then it is easy 
to show that Eq. is valid for low-yield donors as well 
upon replacing 7^ with the free space fluorescence rate 
Id- 



To highlight the role of PERT in the far field RET, 
consider energy transfer from remote donors to an ac- 
ceptor located near the metal surface. In this case, the 
donor's decay rate and spectral function are unaffected 
by metal and RET is dominated by the following process: 
A donor first radiatively excites SP in the metal which 
then nonradiatively transfers its energy to the acceptor. 
The coupling D a d can be derived from Dyson's equation 
for Green's dyadic, 

G(r,r') = G°(r,r') + k 2 lJ dU m G°(r,r m ) • G(r m ,r'), 

where integration is restricted to metal region and e(u) — 
e(w)/eo — 1. For remote donors, using the far field limit 
(kr 3> 1 and kr' <C 1) of the free Green's dyadic^ 

G°(r,r') = ^(6^ u - r M f„), and averaging out over 
donors angular positions and their dipoles' orientations, 
we obtain PERT per donor 
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is the SP enhancement factor for a metal nanostructure 
of general shape. If the acceptor is located at distance r a 
from the center of a spherical NP, we get A = A 1 - cos 2 <f>+ 
A^ sin 2 <b, where 



A = 



Oil 

1 + 2^- 



A» = 



(21) 



are enhancement factors for normal and parallel dipole 
orientations^ oti(uj) is the NP dipole polarizability, and 
cos0 = f a • e a . Eq. (fTSI extends the far field radiative 
RET theory^i 2 ^ to plasmonic systems. In fact, the PERT 
mechanism can dominate RET even in the near field, as 
our numerical calculations below demonstrate. 



III. NUMERICAL RESULTS FOR NEAR FIELD 
ENERGY TRANSFER 

As an example, consider a donor and acceptor near 
spherical Ag NP in water with normal dipole orienta- 
tions (see Fig. 1). The near field matrix Djk is read- 
ily obtained from the Mie's theory Green's dyadic^ 2 - as 
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D 
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u ik 
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FIG. 2: (Color online) (a) Spectral function Eq. (fT6)l and 
molecules' optical bands relative to SP band oti/R 3 (inset) are 
shown together with (b) quenching factor y^/Fd and coupling 
|D a[ j| 2 (inset) using full and nonradiative (NR) models. 



FIG. 3: (Color online) RET vs distance for R = 20 nm Ag 
NP is shown (a) at 9 = n/3 with d a = dd and d a — 2 nm 
(inset) and (b) with d a = 10 nm at 9 — it/ 3 and 8 = n (inset) 
using full, nonradiative (NR), and Forster models. 



are NP-induced radiative and nonradiative terms, a; = 
R 2l+1 iJ+(i~+i) eQ is NP polarizability, Pi(x) is a Legendre 

polynomial, D° d = (1 + sin 2 0/2)/t^ d , f • r d = cos6>, 
and angular momenta up to I = 50 are included. Full 
decay rates are Tj = — (3/2fc 3 ) JjD'jj. We consider, for 
simplicity, a high-yield donor with a broad emission band 
due to the vibrational modes. Molecules optical bands 
are Lorentzians of width 0.05 eV centered at 2.95 eV 
and 3.2 eV with maximal overlap at about SP energy of 
3.08 eV [see inset in Fig. Ufa)], er a (w) was normalized to 
its total J duxj a {uj), and modified a a , fd, and D a d were 
found using Eq. (j22"j) . 

In Fig. [TJ we plot W a d vs the molecule's distance d 
from the R = 30 nm NP surface at 9 = ir/3 with equal 
d a = dd = d and with changing dd = d at fixed da- 
Three models — the full Eq. ([2]), its nonradiative part 
only, and the Gersten-Nitzan mode l 30 ' 31 — are compared 
to Forster's transfer Eq. ([1]). For dd = d a , W a d is about 
three times larger than W^ d and rapidly decays with d, 
while for d/R 1 it is quenched by metal. There is no 
enhancement if only the nonradiative channel is included 
in Eq. ©. In contrast, the Gersten-Nitzan model yields 
much greater enhancement (up to 10 5 ) for d/R <C 1 since 



it includes no quenching effects. However, at fixed d a and 
dd/R > 1, the full Wad is the largest one [see Fig. [TJb)] 
due to the dominant role of the PERT mechanism, as 
discussed above. 

The interplay of different RET contributions is shown 
in Fig. [5] featuring spectral density Eq. (fTo| together 
with quenching factor 72/Td and coupling \D ad \ at fixed 
d. dWad/duJ has a sharp SP peak which disappears if 
only the nonradiative channel is included [see Fig. [Ha)]. 
PERT channel reduces ^/^d due to SP-enhanced radi- 
ation but it strongly enhances D a d [see Fig. [UJb)], the 
net result being RET increase, while in the nonradiative 
channel the enhancement and quenching effects nearly 
cancel out. Weak high-frequency oscillations are due to 
high-; SPs. 

The relative rates of SP-assisted RET and Forster's 
transfer are highly sensitive to the system's geometry. 
RET is quenched if both molecules are close to the NP 
surface [see Fig. 02a)] but it becomes enhanced if donor- 
NP distance increases (inset). For 9 — tt/3 RET is en- 
hanced if dd > R [see Fig. E^b)], but for 9 = tt it is 
strongly enhanced for nearly all d (inset). In fact, NP 
acts as a hub that couples equally well nearby and re- 
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mote molecules with different 6 while Forster's transfer 
drops for large r a( j. For smaller NP sizes, the role of 
PERT becomes less pronounced yet remains dominant 
for larger donor-NP distances. 

IV. CONCLUSIONS 

In summary, a theory of resonance energy transfer 
between energy donors and acceptors near a plasmonic 



structure is presented which maintains a correct energy 
balance between transfer, dissipation, and radiation that 
is essential for interpretation of experimental data. The 
plasmon-enhanced radiative transfer is shown to be the 
dominant mechanism in a wide parameter range. This 
work was supported by the NSF under Grants No. DMR- 
0906945 and No. HRD-0833178, and by the EPSCOR 
program. 
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